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Abstract: We study the problem of addition and subtraction using the Zeckendorf representation of in- 
tegers. We show that both operations can be performed in linear time; in fact they can be performed by 
combinational logic networks with linear size and logarithmic depth. The implications of these results for 
multiplication, division and square-root extraction are also discussed. 



1. Introduction 

Zcckendorf [Z] observed that every integer X > can be represented as the sum of a unique subset of 
the Fibonacci numbers {F n : n > 2} in which no two consecutive Fibonacci numbers appear. That is, we 
may write 

x = y~] x k F k 

k>2 

for a unique sequence Xk such that Xk S {0, 1} and XkXk+i = for all k > 2. This representation is called 
the Zeckendorf representation (or sometimes the Fibonacci representation) of X. Zcckendorf representations 
have applications in coding; they have been proposed for use in self-delimiting codes by Apostolico and 
Fraenkcl [A] and by Fracnkcl and Klein [F2]. They can also be used for run- length- limited binary codes (in 
which neither nor 1 may appear more than twice in succession). In both of these applications, the main step 
of encoding is to convert an integer from its binary representation to its Zeckendorf representation (with 
the reverse conversion occurring at the receiving end). Arithmetic involving integers in their Zeckendorf 
representations also occurs in an algorithm for playing Wythoff's game [W3] due to Silber [S]. 

The main problem we address in this paper is: given the Zcckendorf representations of two integers, X 
and Y , how can we find the Zeckendorf representation of their sum Z = X + Y? Zcckendorf addition, and 
other arithmetic operations in Zeckendorf representation, have been discussed by Ligomenides and Newcomb 
[L], by Freitag and Phillips [F3, F4], and by Fenwick [Fl], but none of these authors explicitly discuss the 
resources required by the methods proposed. Tee [T] reviews these methods and gives explicit bounds, but 
his bounds are extremely weak: his algorithm for Zeckendorf addition of two n-digit numbers runs in time 
0(n 3 ). In Section 2, we shall show that Zeckendorf addition of n-digit numbers can be performed in time 
0(n); in fact, it can be carried out in three linear passes (in alternating directions) over the input sequence. 
This algorithm lends itself to efficient parallel as well as serial implementation; in fact, it can be carried out 
by combinatorial logic networks (acyclic interconnections of gates with bounded fan-in) having size (number 
of gates, used as an estimate of cost) 0(n) and simultaneously having depth (number of gates on the longest 
path from an input to an output, used as an estimate of delay, or parallel execution time) O(logn). Apart 
from the constants implicit in the O-notation, this result is the best possible (because the computation of any 
function (such as the most significant bit of the output) that depends on all the inputs digits requires at least 
linear size and at least logarithmic depth). See Wegener [W2] for an excellent treatment of combinational 
logic networks as a computational model well suited to the discussion of resource bounds. 

In Section 3, we shall extend our results on addition to subtraction (that is, to the addition of signed 
integers, represented by adjoining a positive or negative sign to their Zeckendorf representations) . Our result 
is that signed integers can also be added using three alternating passes over the input, so that all the bounds 
derived in the preceding section continue to hold, albeit with larger constants in the O-terms. 

Finally, in Section 4, we discuss the implications of the foregoing results to the problems of multiplica- 
tion, division (with remainder), and square- root extraction (with remainder). For all three problems, there 
are combinational logic networks treating n-digit numbers with size 0(n 2 ) and depth 0((logn) 2 ). Our 
algorithms for these problems are based on conversion from Zeckendorf to binary representation (for which 
we describe networks of size 0(n 2 ) and depth O(logn)), and from binary to Zeckendorf representation (for 
which we describe networks of size 0{n 2 ) and depth 0((logn) 2 )). 
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2. Addition 

In this section, we shall show how finite automata making three passes over the input can perform 
addition of Zcckendorf representations. We assume that the problem is presented as a sequence of O's, l's 
and 2's, obtained from the Zeckendorf representations of the numbers to be added by adding the digits in 
each position independently. This sequence will in general not be the Zeckendorf representation of the sum, 
since (1) it may have consecutive l's, and (2) it may contain 2's. But we note that any 2's in this sequence 
are both immediately preceded and immediately followed by O's (because the l's in the numbers to be added 
that produced these 2's could not have other l's in consecutive positions). 

We shall compute the Zeckendorf representation of the sum from this sequence in two stages. In the 
first stage we shall eliminate the 2's, obtaining an intermediate sequence that contains only O's and l's, but 
that may contain consecutive l's. In the second stage we shall convert this intermediate sequence into the 
Zeckendorf representation of the sum. 

The first stage will be performed in one left-to-right pass over the input sequence. (We assume that 
the most-significant digits are on the left and the least-significant are on the right, as in conventional radix 
representations. We also assume that the sequence begins with a 0, which may be appended if necessary with 
out changing the value represented.) We describe the actions taken by the algorithm in terms of a "moving 
window" four positions wide. The window begins at the leftmost four positions. At each step it may change 
the values of the digits at positions within the window, in accordance with rules presented below. After 
each step the window is shifted one position to the right. When the window has reached the rightmost four 
positions, and any changes applicable to those positions have been made, a final "cleanup" operation will be 
performed. 

The rules describing the changes to be made at each window position are as follows. 

020a; i — y 100a;' 
030a; >->■ 110a;' 
021a; ^ 110a; 
012a; ^ lOlx 

Here x denotes one of the symbols 0, 1 or 2, and x' denotes its successor: 1, 2 or 3. These rules employ 
the symbol 3 in addition to 0, 1 and 2, but at the end of the first stage all 3's, as well as all 2's, will have 
been eliminated. If none of these rules are applicable, the symbols within the window are not changed. 

The soundness of these rules (that is, the fact that they leave the value represented unchanged) follows 
immediately from the recurrence for the Fibonacci numbers. It remains to prove their effectiveness (that is, 
that all 2's and 3's are eventually eliminated). This effectiveness is a consequence of the following two facts. 

Fact 1: Whenever a 3 is created (which happens when a 2 is incremented in the fourth position of the 
window), it is both preceded by a and followed by a when it reaches the second position of the window, 
and is thus eliminated at that time. 

Fact 2: Every 2 (whether present in the input sequence or created by incrementing a 1 in the fourth position) 
is either (1) preceded by a and followed by a or a 1 when it reaches the second position, and is thus 
eliminated at that time, or (2) is preceded by the sequence 01 when it reaches the third position, and is thus 
eliminated at that time. 
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Proof of Fact 1: A 3 is created by incrementing a 2 in the fourth position. Any symbol that is incremented 
must be preceded by a 0. A 2 that is incremented must have been present in the input, where it must have 
been followed by a 0, and this cannot be incremented because it is not preceded by another 0. □ 

Proof of Fact 2: Consider first a 2 that is created by incrementing a 1 in the fourth position. Any symbol 
that is incremented must be preceded by a 0. A 1 that is incremented must have been present in the input, 
where it must have been followed by either a or a 1, and this or 1 cannot be incremented because it is 
not preceded by another 0. Thus created 2's are eliminated when they reach the second position. 

Consider then a 2 that is present in the input, where it is both preceded and followed by 0. If the 
preceding is not incremented when it is in the fourth position, then the reasoning above applies. If the 
preceding is incremented, then it must have been preceded by another 0, so that the sequence 01 then 
precedes the 2. The that initially followed the 2 cannot be incremented because it is not preceded by 
another 0. Thus this 2 will be eliminated when it reached the third position. □ 

After any applicable changes have been made with the window in its rightmost position, there may still 
be a 2 or 3 in the third or fourth position of the window. These may be cleaned up as follows. If there is 
a 3 in the third position, then it must be preceded and followed by 0's. Then 030 can be changed to 111 
without changing the value represented. If there is a 2 in the third position it must be preceded either by a 
or by the subsequence 01, and it must be followed by a 0. Then 020 can be changed to 101, or 0120 can be 
changed to 1010 without changing the value represented. If there is a 3 in the fourth position, then it must 
be preceded by a 0. Then 03 can be changed to 11 without changing the value represented. If there is a 2 
in the fourth position it must be preceded either by a or by the subsequence 01. Then 02 can be changed 
to 10, or 012 can be changed to 101 without changing the value represented. After this cleanup operation, 
the resulting intermediate sequence consists entirely of 0's and l's. 

The second stage will be performed in two passes over the intermediate sequence, the first from right 
to left, and the second from left to right. The first, right-to-left pass will use a window of width three, and 
will make changes according to the following single rule. 

011 i Y 100 

The soundness of this rule is again follows immediately from the recurrence for the Fibonacci numbers. After 
this right-to-left pass, the resulting sequence contains no occurrence of the subsequence 1011. Suppose, to 
obtain a contradiction, that the resulting sequence does contain an occurrence of the subsequence 1011. 
Focus on one such occurrence, and suppose that at step s the window was positioned at the final three 
positions of this occurrence. At the outset of step s, the final two positions of the occurrence must have been 
1, because at and after step s they could only be changed from l's to 0s, not from to l's. Since they were 
not changed to 0's at step s, the second position of the occurrence must have been 1 at step s, and thus 
must have subsequently been changed to 0. But there are only two steps that could have made this change: 
steps s + 1 and s + 2. If step s + 1 had made the change, it would also have put a in the third position 
of the occurrence, which would have remained forever after, contrary to what we see now. And if step s + 2 
had made the change, it would also have put a in the first position of the occurrence, which would have 
remained forever after, contrary to what we see now. Thus we have reached a contradiction, proving there 
is no occurrence of the subsequence 1011 after the right-to-left pass. 

The second pass of the second stage is a left-to-right pass using the same width-three window and the 
same rule. It is easy to see that during this pass every pair of consecutive l's is eliminated, no new pairs of 
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consecutive l's are created, and no occurrence of 1011 is created. Thus at the conclusion of stage two, the 
resulting sequence is the Zeckendorf representation of the sum. 

Since the operations performed in each pass can be carried out by a finite automaton making a single 
pass over the input, they can also be carried out by combinational logic networks having size O(n) and depth 
O(logn), where n is the length of the input sequence. 

3. Signed Numbers 

In this section, we shall extend the results of the preceding section to cover signed integers. We shall 
assume that a sign bit (0 indicates non-negative, 1 indicates non-positive) is appended to the Zeckendorf 
representation. Thus subtraction can be performed by flipping the sign of the number to be subtracted, then 
performing addition of signed numbers. 

Consider the addition of two signed numbers. If the numbers have the same sign, one adds the mag- 
nitudes of the numbers and appends the common sign to the sum. This can be accomplished by trivial 
extension of the algorithm presented for unsigned addition. If, on the other hand, the numbers have oppo- 
site signs, we must subtract the smaller magnitude from the larger magnitude, and append the sign of the 
number with larger magnitude. There are two problems here: one is to determine which number has the 
greater magnitude, and the other is to determine the difference between the magnitudes. 

We shall describe the solution to the second problem first, assuming that we know that the positive 
number has greater magnitude. After presenting this solution, we shall indicate how to modify it to solve 
the first problem as well. 

We shall assume that the input is represented as a sequence of 0's, +l's and —l's, obtained from the 
Zeckendorf representations of the two numbers by combining the digits in each position separately, canceling 
+l's and —l's that occur in the same position to 0's. For typographical convenience, we shall write +l's 
and —l's as l's and l's, respectively. 

Our algorithm will begin with a preliminary left-to-right pass over the input sequence. The output 
from this pass will be a sequence of 0's l's and 2's that can be used as input to the algorithm described in 
the preceding section to compute the difference between the magnitudes. (Since the latter algorithm also 
begins with a left-to-right-pass, these two passes can be combined into a single pass with a wider window, 
so the algorithm for signed addition can also be performed in three alternating passes.) The preliminary 
left-to-right pass will use a window of width three and apply the following rules. 
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The strategy of these rules is clear. They keep a symbol with positive sign in the window at all times, and 
use it to cancel symbols with negative sign whenever necessary. They may need to introduce 2's, but any 
2's in the output sequence are both preceded and followed by O's. Thus the output sequence satisfies the 
conditions required for the application of the addition algorithm in the preceding section, which can then be 
used to determine the difference between the magnitudes. 

We must now address the problem of determining which of the signed numbers to be added has the 
greater magnitude. To do this, we look for the first occurrence of a 1 or 1 as we scan the input sequence from 
left-to-right; the sign of this occurrence will be the sign of the sum. (If all the symbols of the input sequence 
are 0, the sum is zero.) One way to implement this process is to extend the window of the preliminary pass 
on the right to a fourth position, which will then be the first to see the leading 1 or 1. (We assume that 
three O's are prefixed to the input sequence so that the fourth position will be visited by every symbol of the 
original input.) If it detects a leading 1 (rather than 1), it can then flip the sign of that and every following 
symbol, and the result of the pass will then be the difference between the magnitudes, to which the sign of 
the leading symbol can be appended to give the final representation of the sum. 

4. Multiplication, Division and Square-Root Extraction 

In this section, we shall show how our results on Zeckendorf addition can be used to give efficient 
algorithms for more complicated arithmetic operations. We begin by observing that if we use our addi- 
tion algorithm in Fenwick's [Fl] algorithm for multiplication, we obtain combinational logic networks for 
Zeckendorf multiplication with size 0(n 2 ). Unfortunately, these networks have depth O(nlogn) (because n 
additions, each requiring depth (logn), are performed serially). While we know of no improvement to this 
size bound, we can improve the depth bound considerably, to 0((logn) 2 ). 

The strategy for the remaining results in this section is based on what Freitag and Phillips [F4] call 
"cut [ting] the Gordian knot". We convert the inputs from Zeckendorf representation to binary, perform the 
operation using a known efficient binary algorithm, and convert the result back to Zeckendorf representation. 
To implement this strategy, we must look at algorithms for converting between Zeckendorf and binary 
representations, and review some known results on binary arithmetic. 

Our first observation is that a an n-bit number in binary can be converted to Zeckendorf representation 
by a combinational logic network of size 0(n 2 ) and depth 0((logn) 2 ). To see this, we observe that we 
can construct a tree of Zeckendorf adders that computes the sum of n n-digit numbers with size 0(n 2 ) and 
depth 0((logn) 2 ). (The tree contains 0{n) adders, each of size 0(n) and depth O(logn), each constructed 
according to the result presented in Section 2. There are O(logn) adders on each path from an input to an 
output. The tree thus has size 0(n 2 ) and depth 0((logn) 2 ) in total.) For each bit of the binary input, we 
feed the Zeckendorf representation of the appropriate power of two into one of the adders if that bit is a 1, 
and zero into that adder if that bit is a 0. These additional networks have size 0(n 2 ) in total and depth 
O(l), so the bounds for the tree apply to the entire network. 

To convert from Zeckendorf to binary is even easier. So called "carry-save" adders allow the tree of 
binary adders to have size 0{n 2 ) and depth only O(logn) (see Wallace [Wl]). Thus the size and depth 
required for conversion are dominated by those required in the binary-to-Zeckendorf direction. 

For multiplication, division and square-root extraction, there are binary algorithms whose size and depth 
requirements are also dominated by those for binary-to-Zeckendorf conversion. (Wallace [Wl] gives networks 
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for multiplication of size 0(n 2 ) and depth O(logn) (again using trees of carry-save adders). He also gives 
networks of size 0(n 2 ) and depth 0((logn) 2 ) for division and square-root extraction.) 

There are much more efficient (in terms of both size and depth) networks known for binary multiplication, 
division and square-root extraction (see Pippenger [P] for a survey). It is a challenging open problem to try 
to match their efficiency with networks for the corresponding operations in Zeckendorf representation, either 
by improving the conversion process that is the current bottleneck, or by finding more efficient algorithms 
for these operation directly in Zeckendorf representation. 
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